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Mica-based  materials,  either  laminated  papers  or  layered  single  crystals,  are  popular  in  solid  oxide  fuel 
cell  applications  as  sealing  components.  Their  interface  and  bulk  leak  paths  can  be  considered  as  slits  with 
various  heights,  lengths  and  widths.  A  hydromechanics  model  was  established  based  on  the  geometric 
assumptions  to  predict  the  influence  of  slit  geometry  and  pressure  difference  of  the  seal  on  H2  leak  rate. 
The  dependence  of  leak  rate  on  slit  geometry  and  the  pressure  difference  as  well  as  the  effectiveness  of 
compressive  loading  and  compliant  interface  layer  were  discussed  accordingly.  The  model’s  applicability 
is  supported  by  reported  experiments. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Among  the  technical  challenges  in  developing  planar  solid  oxide 
fuel  cells  (SOFCs),  the  sealing  material  has  been  regarded  as  one 
of  the  most  significant  issues.  The  seals  are  required  to  be  non- 
conductive,  and  thermally,  chemically  and  mechanically  stable  in 
both  oxidizing  and  reducing  environments,  and  prevent  fuel  and 
oxidant  gas  leak  within  the  stack  [1-3]. 

The  commonly  used  seals  for  the  planar  SOFCs  have  been  glasses 
or  ceramic-glasses  [2].  The  advantage  of  glass  and  ceramic-glass 
seals  is  that  their  compositions  can  be  tailored  to  optimize  the 
required  physical  properties,  such  as  the  coefficient  of  thermal 
expansion  (CTE);  however,  they  tend  to  transform  in  phases  and 
react  with  the  cell  component  and  interconnect  materials  under 
SOFC  operating  conditions  in  a  long  run,  due  to  their  intrinsic  ther¬ 
modynamical  instability  [4-6]. 

Compressive  seals  have  been  developed  to  avoid  the  disadvan¬ 
tages  of  the  rigid  glass-based  seals,  with  the  merit  of  flexibility  and 
compressibility;  a  close  match  of  CTE  to  those  of  the  adjacent  mate¬ 
rials  is  not  required,  allowing  the  cells  and  interconnects  to  expand 
and  contract  freely  during  thermal  cycles  and  operation,  and  the 
replacement  of  the  malfunctioning  components  became  possible 
[7].  So  far,  two  kinds  of  compressive  seals  have  been  considered, 
i.e.,  the  deformable  metallic  seals  and  the  mica-based  seals.  The 
deformable  metallic  seals  include  ductile  silver  [8,9]  and  corru¬ 
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gated  or  C-shaped  superalloy  gaskets  [10].  Its  application  is  limited 
by  its  high  electronic  conductivity  [9,10].  The  plain  mica  seals  and 
hybrid  mica-based  seals  have  been  systematically  investigated  and 
their  sealing  properties  were  reported  [11-20].  The  plain  mica 
has  exhibited  very  poor  sealing  characteristics,  and  in  contrast, 
the  hybrid  mica-based  seals  have  demonstrated  excellent  hermetic 
property  and  are  being  regarded  as  the  potential  candidate  for  inter¬ 
mediate  temperature  planar  SOFCs.  It  has  been  identified  [7,17,19] 
that  there  are  two  kinds  of  possible  leak  paths  for  mica-based  seals: 
one  is  through  the  body  of  the  mica  (bulk  leak)  and  the  other  is 
through  the  sealing  interfaces  between  the  mica  and  the  adjacent 
SOFC  components  (interface  leak).  Such  interfaces  are  considered 
to  be  the  major  leak  paths  of  the  conventional  mica  seals.  The  pur¬ 
pose  of  the  present  study  is  to  develop  a  hydromechanics  model 
to  predict  the  leak  rate;  and  accordingly,  the  factors  that  control 
the  leak  rate  can  be  identified  and  analyzed,  and  appropriate  seal¬ 
ing  design  procedures  can  be  suggested  to  improve  the  hermetic 
property  of  the  mica-based  seals. 

2.  Modeling 

In  planar  SOFCs,  H2  is  usually  used  as  the  fuel,  and  the  sealing 
requirement  is  higher  than  other  gases  due  to  the  small  molecule 
size  of  H2.  In  the  present  study,  H2  is  chosen  as  the  subject  of  investi¬ 
gation,  ideal  gas  assumption  is  applied  owing  to  its  high  rarefaction. 

It  has  been  observed  by  a  scanning  electron  microscope  (SEM), 
that  the  leak  paths  of  mica  seals,  either  muscovite  single  crystal 
and  paper  [11]  or  phlogopite  paper  [16],  are  mainly  composed  of 
tiny  openings  located  either  at  the  sealing  interfaces  or  in  between 
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Cel I 

Interconnect 

H -  L  - H 

Fig.  1.  Schematic  illustration  of  sealing  configuration  in  a  planar  SOFC  stack. 
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Fig.  2.  Flow  field  set  up  in  Cartesian  coordinates. 


mica  layers.  Both  leak  paths  can  be  idealized  as  the  slits  with  spe¬ 
cific  geometry,  as  schematically  shown  in  Fig.  1.  L  is  the  slit  length, 
which  is  more  or  less  related  to  the  sealing  width.  The  slit  may  not 
be  smooth,  and  the  surface  roughness  is  likely  to  be  in  the  same 
order  of  magnitude  as  that  of  the  molecular  free  path  of  H2.  H2 
leaks  through  the  slits  once  the  pressure  inside  of  the  stack  (inner 
pressure)  Pj  is  bigger  than  the  pressure  of  outside  of  the  stack  (exter¬ 
nal  pressure)  P0  that  equals  to  the  atmospheric  pressure  (see  Fig.  1 ). 
Since  the  flow  is  in  such  a  confined  channel,  the  flow  characteristics 
of  the  high-temperature  gas  should  be  considered. 

In  hydromechanics,  the  non-dimensional  parameter  Kn ,  the 
extent  of  gas  rarefaction,  is  defined  as  Kn  =  A/D ,  where  A  is  the  mean 
free  path  of  gas  molecules  and  D  is  the  flow  dimension  that  herein 
represents  the  slit  height  h.  The  value  of  Kn  determines  the  nature 
of  the  flow.  For  the  continuum  flow,  Kn  is  less  than  0.1 ;  the  slip  flow 
has  a  value  of  Kn  in  the  range  of  0.01  -0.1 ;  for  the  transition  flow, 
Kn  is  in  between  0.1  and  ~10;  and  for  the  free-molecular  flow,  Kn 
is  larger  than  10. 

The  mean  free  path  of  gas  molecules  A  can  be  expressed  as 


condition  [21]  is 

_  h2  dP  Kn 
5  ~  2/z  dx  TTKn 


(2) 


where  /z  is  the  viscosity  coefficient.  Considering  the  leak  flow  of  H2 
in  a  planar  SOFC,  further  pre-conditions  are  given  as  follows: 


(1)  All  the  flows  are  isothermal  and  steady. 

(2)  Temperature  gradient  is  neglected  inside  the  leak  flow  of  H2. 

(3)  The  inertia  term  is  ignored,  which  is  much  smaller  compared 
to  the  diffusion  term  in  the  momentum  equations,  due  to  the 
Reynolds  number  of  the  flow  is  low  in  channels  with  a  large 
length/height  ratio  (L/h  » 1 ). 

(4)  The  rarefaction  coefficient  Cr(Kn)  is  given  as 

Cr{I<n)  =  1  +al<n  (3) 

where  a  =  1.358(2/7r)tan-1  (2.2  Kn0-5 )  [21  ]. 


A  = 


kT 

V2nd2P 


(1) 


where  k  is  Boltzmann  constant,  T is  absolute  temperature,  d  is  effec¬ 
tive  diameter  of  gas  molecular,  and  P  is  gas  pressure.  The  average 
molecular  free  path  of  H2  gas  is  0.63  |xm  in  a  planar  SOFC  operated 
at  800  °C.  If  a  slit  height  (ft)  range  of  0.2-20  p,m  is  assumed  accord¬ 
ing  to  the  observation,  the  pressure-driven  leak  flow  of  H2  gas  falls 
into  the  regime  of  the  slip  flow  or  the  transition  flow. 

In  the  slip  flow  regime  and  the  transition  flow,  the  usually 
assumed  no-slip  boundary  condition  is  no  longer  suitable  and  the 
Knudsen  layer  with  a  thickness  of  approximately  a  molecular  free 
path  plays  a  decisive  role;  in  addition,  the  collisions  between  the 
gas  molecule  and  the  slit  wall  become  more  and  more  important 
compared  to  the  collisions  among  gas  molecules;  and  the  kinematic 
viscosity  needs  to  be  modified  for  the  enhanced  rarefaction  effect. 
Therefore  a  unified  flow  model  for  both  the  slip  and  transition  flow 
regimes  is  to  be  developed,  in  which  the  flow  in  the  slit  is  assumed 
to  be  divided  as  the  Knudsen  layer  and  the  middle  layer,  the  flow 
transition  from  the  Knudsen  layer  to  the  middle  layer  is  continu¬ 
ous,  and  the  enhanced  rarefaction  effect  is  taken  into  account  by 
introducing  the  rarefaction  coefficient  Cr(Kn),  a  function  of  Kn. 

The  flow  in  the  Knudsen  layer  does  not  follow  the  Navier-Stokes 
equation;  a  specific  solution  from  the  Boltzmann  equation  is 
needed.  Flowever,  it  is  intricate  to  solve  the  Boltzmann  equa¬ 
tion,  only  some  semi-analytic  solutions  or  numerical  solutions 
can  be  derived  for  special  cases.  Therefore,  it  is  necessary  to  sim¬ 
plify  the  solving  process  without  changing  the  dependence  of  the 
velocity  slip  boundary  condition  on  Kn.  Usually,  for  a  relatively 
accurate  prediction  of  gas  flow  rate  (leak  rate),  second  order  pre¬ 
cision  velocity  slip  boundary  condition  is  enough.  So  the  general 
velocity  slip  boundary  condition  was  adopted.  For  the  pressure- 
driven  incompressible  flow  in  the  slits,  the  velocity  slip  boundary 


Based  on  the  above  assumptions,  the  velocity  distribution  in  the 
slit  can  be  acquired,  and  finally  the  mass  flow  rate  per  flow  width 
is  derived  by  the  integral  of  the  flow  velocity. 

It  is  assumed  that  the  leak  flow  is  a  steady  state  flow,  indepen¬ 
dent  of  time.  The  streamline  direction  was  set  along  the  x-axis  of 
the  Cartesian  coordinates  as  shown  in  Fig.  2,  consequently,  the  flow 
velocity  has  only  a  component  u  in  the  x-direction;  in  the  y  and 
z  directions,  the  components  are  zero,  that  is,  v  =  w  =  0,  and  the 
continuity  equation  for  incompressible  fluids  becomes 


Since  u  is  constant  in  the  z-direction,  i.e.,  9u/9z  =  0,  hence,  u  is  only 
a  function  of  y,  u  =  u{y). 

If  the  viscosity  coefficient  fi  of  H2  is  constant,  Navier-Stokes 
equations  can  be  written  as 


dP 

=  PSx~Tx+p 


f  d2U  d2U  d2u\ 

^9x2  +  dz2  J 


dP 

=  pgy-^y+/i 


f  d2v  d2v  92v\ 

^9x2  +p+  a?  J 


dP 

■  PSz~  dz+P 


( 92w  92w  92w\ 

(  9x2  +  dy2  +  dz2  J 


(4) 

(5) 

(6) 


Ignoring  the  gravity  force  of  H2  ,gx=gy=gz  =  0,  the  above  equations 
are  simplified  as 


„  dP  d2u 
0  —  — "a — f  P  9 
9x  9y2 


(7) 


(8) 
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o  =  4 


(9) 


P  is  only  the  function  of  x.  Since  the  pressure  difference  between  P; 
and  P0  is  small,  dP/dx  can  be  considered  as  a  constant,  and  u  can  be 
obtained  by  integral  of  Eq.  (7). 


(10) 


u~2 )4y2+Ciy  +  C2 

Applying  the  boundary  conditions  u{ 0)  =  u{h)  =  Us ,  the  velocity  dis¬ 
tribution  in  the  streamline  direction  is  derived  as 


u{y)=!]fM(y^_y  JO i_ 
{y)  2/xd  x\h2  h  \+l<n 


(11) 


The  volume  flow  rate  in  an  unit  cross-section  can  be  obtained 
accordingly  as 

rh 


*-/ 


...  ..  h3  dP  f.  6 Kn 

u(y)dy=-uji^r  +  TTm 


(12) 


Considering  the  effect  of  enhanced  rarefaction  and 
dP/dx  =  -(Pj  —  P0)/L  (L  is  the  depth  of  the  slit),  Eq.  (14)  is  modified 
as 


Q  = 


h3  Pi-Po  r 


12  ii 


1  + 


6Kn 


1  +/Cn 


Cr(I<n) 


(13) 


And  therefore,  the  associated  mass  flow  rate  of  a  single  slit  with  a 
unit  width  can  be  expressed  as 

2P0  h3  Pi-Po 


M  =  1  pQ  = 


RT  12  ii 
P0(Pi-P0)h3 


1  + 


6Kn 


+12 


12  iiRTL 
a  -  1 


1  +  Kn 
[{n  +  1 )  +  2(6  +  a)I<n 
TJ  +  Kn 


Cr(/Cn) 


77-1 


Kn1  In 


f  n  +  Kn\~ 

V  1  +I<n  ). 


(14) 


where  p  =  2 P0/PT  and  Y[  -  Pi /P0 - 

In  reality,  there  are  numerous  (nz)  leak  paths  similar  to  Fig.  1  in 
a  mica-based  seal  with  various  slit  height  hz  and  width  bx  (i  =  1, 2, 3, 
...,/<-  1,  /<),  the  total  leak  mass  flow  rate  is  the  weighted  summation 
of  Mi  in  terms  of  slit  number  and  width. 


Mi 


total 


(15) 


1=1 


3.  Discussions 


The  leak  mass  flow  rate  expressed  by  Eq.  (14)  for  one  slit  shows 
that  M  is  a  function  of  slit  height  h  nearly  up  to  the  third  power. 
Fig.  3  demonstrates  the  dependence  of  H2  leak  rate  at  1073  K  on  h 
with  Pj  as  1.048  x  105  Pa,  /jl  as  1.92  x  10-7  Pa  S,  and  the  slit  length  L 
as  1  cm.  The  H2  leak  rate  decreases  sharply  with  the  decreasing  of 
h:  M  is  4.51  x  lCT8,  8.19  x  lCT6  and  1.56  x  lO^gs-1  cm'1  when  h 
is  0.5, 5  and  15  |jim,  respectively,  30  times  increase  in  h  results  in  ca. 
3500  times  increase  in  leak  rate.  This  reflects  the  strong  impact  of 
Knudsen  flow  and  the  associated  rarefaction  effect  on  the  mass  leak 
rate,  the  smaller  the  slit  height,  the  more  significant  the  influence 
of  the  Knudsen  flow. 

According  to  Eq.  ( 14),  M  is  proportional  to  1  /L  when  other  param¬ 
eters  remain  unchanged.  The  slit  length  is  somewhat  related  to  the 
sealing  width.  If  a  single  slit  penetrates  the  entire  sealing  width,  L 
is  equal  to  the  sealing  width;  however,  in  a  real  mica-based  seal,  a 
single  leak  slit  usually  can  only  penetrate  a  portion  of  the  sealing 
width,  the  actual  leak  path  is  proportional  to,  but  much  longer  than 
the  sealing  width,  and  the  leak  rate  will  be  dramatically  decreased 
by  increasing  the  sealing  width.  In  a  specific  stack  design,  increasing 


h,  |jm 

Fig.  3.  Model  prediction  of  H2  mass  leak  rate  as  a  function  of  slit  height. 


sealing  width  may  cost  the  active  area  and  power  output,  therefore, 
the  effectiveness  of  lowering  the  leak  rate  by  increasing  leak  path 
L  is  limited  in  real  applications. 

For  a  specific  leak  slit  with  a  width  of  bz,  the  contribution  of  the 
slit  width  to  the  mass  leak  rate  is  demonstrated  in  Eq.  (15)  as  a 
factor  of  the  unit  width  leak  rate  M.  From  Eq.  (15),  it  also  can  be 
seen  that  the  total  mass  leak  rate  is  proportional  to  the  leak  slit  size 
distribution  ( n\ ,  n2,  n3, . . .,  nz, . . .,  nk )  and  the  total  leak  slit  number 
N  (N  =  ni)-  With  the  same  total  slit  number  N,  the  total  mass 
leak  rate  may  be  different,  depending  on  the  distribution  (ni,  n2, 
n3,  . . .,  ri;,  . . .,  nk).  On  the  other  hand,  a  smaller  N  may  result  in  a 
bigger  leak  rate  than  a  bigger  N,  if  the  slit  size  distribution  (rii,  n2, 
n3,  . . .,  ?+,  . . .,  nk)  represents  more  large  leak  slits,  especially  with 
bigger  h. 

In  order  to  drive  the  gas  flow,  the  gas  pressure  inside  the  SOFC 
stacks  is  usually  1.4-3.51<Pa  (0.2-0.5psi)  higher  than  the  outside 
atmosphere  pressure,  and  may  be  increased  to  1.4  x  104  Pa  (2  psi) 
with  the  fluctuation  of  source  gas  pressure.  Fig.  4  shows  the  H2 
leak  rate  with  different  Px  and  constant  h  (1  p,m)  and  L  (1  cm). 
The  H2  mass  leak  rate  decreases  linearly  with  (Pj  -P0),  which  is 
consistent  with  the  previous  study  [11  ]  where  the  leak  rates  of  phl- 
ogopite  paper,  muscovite  paper  and  muscovite  single  crystal  were 
measured  against  the  gauge  pressure. 

Based  on  the  analysis  above,  two  suggestions  can  be  made  for 
a  tighter  sealing:  (1)  increasing  the  compressive  load;  and  (2) 
applying  a  compliant  layer,  such  as  soft  metal  or  glass,  at  the 
sealing  surfaces.  With  higher  compressive  stress,  the  slit  height 


Fig.  4.  Model  prediction  of  H2  mass  leak  rate  as  a  function  of  pressure  difference 
between  inside  and  outside  of  SOFC  stack. 
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h  and  width  b  tends  to  be  smaller  and  some  tiny  leak  slits  may 
be  completely  closed,  resulting  in  a  lower  leak  rate.  Jeffrey’s  work 

[2]  showed  that  the  leak  rate  for  all  kinds  of  the  mica  seals,  no 
matter  mica  papers  or  single  crystal  layers,  decreased  with  increas¬ 
ing  applied  compressive  stress,  which  is  supportive  to  the  model 
prediction.  For  muscovite  single  crystals,  increasing  compressive 
loading  is  more  effective  due  to  its  smoother  slit  surfaces  and  readi¬ 
ness  for  slit  closure;  a  leak  rate  of  1/10  to  ~l/5  of  that  from  the 
mica  papers  was  achieved  at  the  same  compressive  pressure  [11]. 
The  robustness  of  mica-based  sealing  can  be  further  improved  by 
using  a  silver  or  specific  glass  at  the  sealing  surface,  which  can 
eliminate  the  major  leak  paths  in  the  interface;  or  by  infiltration 
of  a  second  phase,  such  as  glass,  to  seal  the  interface  and  bulk 
leaks,  a  leak  rate  down  to  1.3  x  10-4  seem  cm-1  can  be  resulted 
[17]. 

4.  Conclusions 

According  to  the  model  study  and  experiment  comparison,  the 
following  conclusions  can  be  made; 

(1)  The  leak  paths  in  mica-based  seals  can  be  simulated  as  slits 
with  various  heights  h,  lengths  L  and  widths  b.  The  leak  rate  is 
a  function  of  the  slit  height  h  nearly  up  to  the  third  power ,  and 
is  proportional  to  the  slit  width  and  the  reciprocal  of  the  slit 
length. 

(2)  The  total  leak  rate  increases  with  the  total  slit  number  of  the 

seal  N,  and  is  strongly  related  to  the  slit  size  distribution  (rii, 
n2,n3 . n; . nk). 

(3 )  Pressure  difference  between  inside  and  outside  of  the  cell  drives 
H2  leaking  through  the  seal,  and  the  leak  rate  is  proportional  to 
the  pressure  difference. 


(4)  Increasing  the  compressive  loading  and  applying  compliant 
interlayer  improve  the  robustness  of  the  seal  by  reducing  the 
leak  slit  size  and  seal  the  leak  paths  at  the  sealing  interfaces. 
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